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Context. Convection in stars excites resonant acoustic waves. The frequencies of these oscillations depend on the sound speed inside the star, 
which in turn depends on density, temperature, gas motion, and other properties of the stellar interior Therefore, analysis of the oscillations 
provides an unrivaled method to probe the internal structure of a star. 

Aims. Solar-like oscillations in the primary of the visual binary 70 Ophiuchi are investigated. 

Methods. 70 Ophiuchi A was observed with the Harps spectrograph mounted on the 3.6-m telescope at the ESO La Silla Observatory (Chile) 
during 6 nights in July 2004 allowing us to collect 1758 radial velocity measurements with a standard deviation of about 1.39ms"'. 
Results. The power spectrum of the high precision velocity time series clearly presents several identifiable peaks between 3 and 6mHz 
showing regularity with a large spacing of Av = 161.7 ± 0.3 fjHz. Fourteen individual modes were identified with amplitudes in the range 1 1 to 
14 cm s"'. 
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1. Introduction 

Analysis of the oscillation spectrum provides an unrivaled 
method to probe the stellar internal structure, since the frequen- 
cies of these oscillations depend on the sound speed inside the 
star, which in turn depends on density, temperature, gas motion, 
and other properties of the stellar interior. The prime example 
of such a probe is the Sun: the five-minute oscillations have 
provided a wealth of information about the solar interior. The 
Sun oscillates simultaneously in many modes with periods of 
about 5 minutes and Doppler amplitudes of about 23 cm s"' for 
the strongest modes. Observation of thousands of oscillation 
frequencies and comparing them with theoretical calculations 
has led to significant revisions of solar models. 

These results stimulated various attempts to detect a sim- 
ilar signal on other solar-like stars and the stabilized spectro- 
graphs have achieved the accuracy needed for solar-like os- 
cillation detection by means of radial velocity measurements 
(Carrier et al. I2003> . However, a major asteroseismological 
difficulty is the confrontation of observations and theoretical 
models (Eggenberger et al. I2004bl ^005 1. The observational 
measurements available for an isolated star - such as the effec- 
tive temperature, the metallicity, the luminosity, and eventually 
the interferometric radius, combined with oscillation frequen- 
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cies - provide strong constraints to the global parameters of the 
star but are often not sufficient to unambiguously determine 
an adequate model and to really test the physics of the mod- 
els. The additional constraints imposed by the binary nature, 
namely the same age and initial chemical composition, are ex- 
tremely valuable for accurately determining the properties of 
a binary system. Moreover, in the case of binaries, the masses 
of both components are accurately known by combining visual 
and spectroscopic orbits. A primary target for the search for 
p-mode oscillations in such a system is the bright KO dwarf 
70 Ophiuchi A (HD 165341A), in addition to the a Cen system 
(Bouchy & Carrier l2002l Carrier & Bourban l2003l Bedding et 
al. "2004 and Kjeldsen et al. 2005 i and Procyon (Eggenberger 
et al. ■2004a. Martic et al.2Q04.l- 

In this paper, we report Doppler observations of 
70 Ophiuchi A made with the Harps spectrograph resulting in 
the detection of p-mode oscillations. The observations and data 
reduction are presented in Sect. |2j the acoustic spectrum anal- 
ysis and the mode identification in Sect. |3] and the conclusion 
is given in Sect.|3 



2. Observations and data reduction 

With a V-magnitude of 4. 1 , 70 Oph A is the main component of 
a spectroscopic visual binary composed of a KO and a K5 dwarf 
with an orbital period of 88.38 yr (Pourbaix 2000 1. It was ob- 
served over 6 nights in July 2004 with the Harps spectrograph 
(Pepe et al. I2002> mounted on the 3.6-m telescope at La Silla 
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Fig. 1. Radial velocity measurements of 70 Oph A. Top: A 
constant value of -5.71kms"' was removed from the origi- 
nal data. Bottom: An 8-order polynomial fit was subtracted 
(dashed line on the top window). 

Table 1. Distribution and dispersion of Doppler measurements. 



Date 


No spectra 


No hours 


o" (m s ' ) 


2004/07/05 


81 


3.58 


1.54 


2004/07/06 


290 


7.78 


1.76 


2004/07/07 


350 


7.80 


1.39 


2004/07/08 


390 


7.95 


1.52 


2004/07/09 


331 


8.07 


1.00 


2004/07/10 


316 


6.45 


1.11 



Observatory (ESO, Chile). We took sequences of 40-60 s expo- 
sures, depending on the airmass and the extinction, with a dead 
time of 30 s in-between. In total, 1758 spectra were collected 
with a typical signal-to-noise ratio (S/N) in the range of 120- 
220 at 530 nm. During the stellar exposures, the spectrum of 
a thorium lamp carried by a second fiber was simultaneously 
recorded in order to monitor the spectrograph's stability. The 
spectra obtained were extracted on-line. The radial velocities 
were computed by weighted cross-correlation with a numerical 
mask constructed from a K5 dwarf spectrum. They were also 
determined by the optimum-weight procedure (Connes 119851 
Carrier et al.'2001't but without significant gain. 

The modulation of the Doppler measurements is due to 
the binarity of the star (see Fig. Q). Thus, an 8-order poly- 
nomial fit was subtracted from the raw data to best eliminate 
the Keplerian movement and variations caused by tempera- 
ture fluctuations of the instrument during a few days due to a 
power failure that happened just before our observation run. 
This high-pass filtering does not affect the results of the p- 
modes analysis: indeed, a 2 or an 8-order polynomial fit lead 
to the same power spectrum above 2 mHz. The rms scatter of 
this "new" time series is 1.39 ms"' (see Table^, the solar-like 
oscillations representing a non-neglecting part of this number. 



3. Power spectrum analysis 

In order to compute the power spectrum of the velocity time 
series, we use the Lomb-Scargle modified algorithm (Lomb 
119761 Scargle'1982V Its time scale gives a formal resolution of 
2.2 juHz. The resulting periodogram, shown in Fig.|2] exhibits a 
series of peaks between 3 and 6 mHz, exactly where the solar- 
like oscillations are expected for this star. Typically for such a 
power spectrum, the noise has two components: 

- At high frequencies it is flat, indicative of the Poisson statis- 
tics of photon noise. The mean white noise level cr^o^ cal- 
culated between 1.5 and 2.5 mHz is 0.0026 m^ s"^, namely 
cTji^p = -\/o"pow * 7r/4 = 4.5cms"' in amplitude (Kjeldsen 
& Bedding I1995> . With 1758 measurements, this high 
frequency noise corresponds to crgv - -\jNcrpov//4- = 
1.07 ms"'. This radial velocity uncertainty is larger than 
the uncertainty due to photon noise, which is estimated at 
0.25-0.4ms This discrepancy could be explained by the 
bad seeing and the technical problems occurring during the 
observation run. However, part of the noise should have a 
stellar origin (like granulation). 

- Towards the lowest frequencies, the power should scale in- 
versely with frequency squared, as expected for instrumen- 
tal instabilities. Moreover, this star shows low frequencies, 
not completely removed by the polynomial fit, due to the 
Keplerian movement. 

We note that a high peak near 3.1 mHz and a smaller one 
nearly twice this value are widely decentered with regard to 
the others (see Fig.|2}. The ESO staff pointed out that a large 
noise peak could occur with a period of about 6 minutes (near 
3 mHz and multiple frequencies at 6 and 9 mHz). This noise 
is not completely understood at present; it could be due to the 
guiding, although this seems unlikely given that the guiding 
correction was applied every 5 seconds, or to the telescope. 
Afterwards, we thus limit our p-modes research between 3.2 
and 6 mHz. 

3.1. Search for a comb-like pattern 

In solar-like stars, p-mode oscillations of low-degree are ex- 
pected to produce a characteristic comb-like structure in the 
power spectrum with mode frequencies v„/ reasonably well ap- 
proximated by the asymptotic relation (Tassoul ll980t : 

v„j ^ Av(n +€)-(((+ \)Do . (1) 

Here, Do, is sensitive to the sound speed near the core and is 
equal to ^6vq2 when the asymptotic relation holds exactly; and 
e is sensitive to the surface layers. The quantum numbers n and 
{ correspond to the radial order and the angular degree of the 
modes, and Av and 6vq2 to the large and small separations. To 
search for periodicity in the power spectrum, an autocorrela- 
tion with a threshold of 80 cm^ s"^ (all values below this limit 
are fixed to zero) is calculated and presented in Fig. |3] Each 
peak of the autocorrelation corresponds to a structure present 
in the power spectrum. The three strong peaks at low frequency 
at about 11.5, 23, and 35/iHz correspond to the daily aUases. 
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Fig. 2. Power spectrum of 70 Oph A. The high peak in gray 
(and twice this value) should be due to a technical problem 
(probably due to a hard point of the telescope). The Nyquist 



frequency has a value of 7.14mHz. The window function is Fig. 3. Autocorrelation ofthe power spectrum of 70 Oph A with 
shown in the inset with daily aliases at 1 1 .57 and 23. 15/iHz. a threshold of 80 cm^ s"^. 



The next dominant peaks appear between 120 and 170/iHz, 
and one of them corresponds to the large spacing. Note that 
the external parameters of 70 Oph A, scaling from the solar 
case (Kjeldsen & Bedding 1995), give a large frequency spac- 
ing between 150 and 170yuHz. Indeed, by refining the orbit of 
the system and measuring the individual magnitudes and effec- 
tive temperatures (for both components), we obtain the follow- 
ing parameters for 70 Oph A: a mass near 0.87 Mq, a luminos- 
ity of 0.52 Lq, an effective temperature of 5300 K, and a solar 
metallicity. Details about these parameters are postponed to the 
theoretical paper (Eggenberger et al. in preparation). Other im- 
portant peaks in the autocorrelation are situated between 320 
and 380jL/Hz, corresponding to twice the large separation. The 
large spacing will then correspond, knowing the domain of one 
time and twice the large spacing, to half of 333.8 + 1 1 .57 yuHz, 
namely the two strong peaks at about 161.5 and 172jt/Hz and 
the small one at 167juHz. To consolidate our solutions we se- 
lected all those peaks between 3.2 and 6 mHz with an amplitude 
greater than 10 cms"' and tried to build echelle diagrams with 
different large separations. Solutions showing well-aligned fre- 
quencies are found for large spacings of about 161.5, 167, and 
172 //Hz. 

3.2. Mode identification 

The method used to extract mode frequencies is described well 
in Carrier et al (2005a 2005b i. The frequencies are extracted 
using an iterative algorithm that identifies the highest peak be- 
tween 3.2 and 6 mHz and subtracts it from the time series. Note 
that because of the stochastic nature of solar-like oscillations, 
a timestring of radial velocities cannot be expected to be per- 
fectly reproduced by a sum of sinusoidal terms. Therefore, us- 
ing an iterative clean algorithm to extract the frequencies can 



add additional peaks with small amplitudes due to the finite 
lifetimes of the modes that we do not know. In the case of 
70 Oph A, we iterated the process until all peaks with an am- 
plitude higher than 2.5 cr in the amplitude spectrum were re- 
moved. Then, cr represents the noise in the amplitude spectrum 
and has a value of 4.5 cms"' (see Sect. E}- This threshold en- 
sures that the selected peaks have only a small chance of being 
due to noise. 

Using the extracted frequencies, we drew echelle diagrams 
with different large spacings. The best echelle diagrams are 
found for the large spacings of 161.7 and 172.2yuHz. The so- 
lution with a large spacing of 167 //Hz, found in the previous 
section, is ruled out by the echelle diagram; indeed, the dis- 
tance between the corresponding { = (or 2) and ^ = 1 modes is 
not compatible with such a large separation. The results of the 
automatic extraction are presented in Table |2] for the two dif- 
ferent possible large separations with the corresponding identi- 
fication. Because of the daily alias of 1 1.57//Hz introduced by 
the mono-site observations, we cannot know a priori whether 
the frequency selected by the algorithm is the right one or an 
alias. We thus considered that the frequencies could be shifted 
by +11.57 //Hz. 

To investigate how many peaks are expected to be due to 
noise, we conducted simulations in which we analyzed noise 
spectra containing no signal. For this purpose, a velocity time 
series was built, using the observational time sampling and ra- 
dial velocities randomly drawn by assuming a Gaussian noise 
(Monte-Carlo simulations). The amplitude spectrum of this se- 
ries was then calculated and peaks with amplitude greater than 
2.5, 3, and 3.5 cr were counted; note that a peak and its aliases 
are only counted once. The whole procedure was repeated 1000 
times to ensure the stability of the results. In this way, we find 
that the number of peaks due to noise with an amplitude larger 
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Table 2. Identification of extracted frequencies. The signs + or 
- mean that the peaks have to be shifted to a higher or lower 
value, respectively, by 1 1.57 yuHz. 



Frequency 


Mode ID (Avi) 


Mode ID (Av2) 


S/N 


L"Hz] 


161.7 yuHz 


172.2 //Hz 




3246.2 


noise 


+ {=1 


3.2 


3350.9 


noise 


{ = 


2.8 


3378.1 


{ = 


noise 


2.6 


3552.1 


-{ = 


noise 


2.6 


3862.1 


( = 


( = 


3.0 


4194.4 


-( = 


+ { = 


3.2 


4258.8 


e = 1 


noise 


2.9 


4343.7 


e = o 


noise 


2.8 


4650.4 


+ e = 2 


-{=1 


2.8 


4719.6 


noise 


{ = 


2.5 


4797.9 


noise 


+ e=i 


3.6 


4993.2 


{ = 


~{=l 


2.5 


5054.4 


+ e=i 


+ e = o 


2.8 


5153.7 


i = 


t = 1 


3.1 


5255.2 


noise 


-{ = 


2.9 


5328.5 


~{ = 


t = 1 


2.6 


5487.3 


-( = 2 


+ e=i 


2.7 


5717.8 


{ = 1 


noise 


2.7 


5763.3 


noise 


{ = 


3.4 


5890.9 


= 1 


noise 


2.6 


5990.0 


noise 


noise 


2.7 




7 noise peaks 


7 noise peaks 





than 2.5 cr is 13.0±4.6 in the range 3.2-6 mHz, for 3cr it is 
1 .9 + 1 .4, and for 3.5 cr the number of peaks due to noise varies 
between and 3 with a mean value of 0.2 and a standard de- 
viation of 0.4. These results are in good agreement with the 
number of identified noise peaks (see Table|2}. 

The star 70 Oph A is very similar to a Cen B (see Carrier & 
Bourban l2003t . the same spectral type and similar large spac- 
ing, which has a mean small spacing of 10 //Hz in the frequency 
range 3-4.5 mHz. Inspecting the results of the mode identifica- 
tion, we note that the value of the small spacing coming from 
the identification with the large separation of 172.2//Hz is sig- 
nificantly different from 10 //Hz. No matter what the identifi- 
cation may be (f = 0, 1 or 2 for both lines in the echelle dia- 
gram), the small separation will be lower than 6.5 //Hz in the 
frequency range 3^.5 mHz. This value of the small spacing, 
which is lower than expected, suggests that this identification is 
less reliable than the one with a large separation of 161.7//Hz. 
Although the solution Av - 172.2//Hz cannot be definitely 
ruled out, we afterwards consider only the one with a large 
spacing of 161.7/(Hz. 

The echelle diagram based on the extracted frequencies 
shows only two lines. It is not possible to disentangle £-2 from 
^ = modes with our time series, mainly due to the poor reso- 
lution (2.2//Hz): daily aliases of ^ = modes merge with i-2 
modes, and vice versa. This suggests that the small spacing is 
close to 10 //Hz. As a result, some extracted { = 2 frequencies 
could be shifted by 1 1 .57 //Hz and be identified as radial modes 
and inversely. 




50 100 150 

modulo AiJ [jtiHz] 

Fig. 4. Echelle diagram of identified modes with a large sepa- 
ration of 161.7//HZ. The modes {^2 (■), i^O (•), and ^=1 (A) 
are represented with a size proportional to their amplitude. 

Due to the ambiguous identification of ^ = 0, 2 modes, an- 
other identification can be found replacing all { = 0,2 by ^ = 1 
modes and all ^ = 1 by ^ = 2 modes. In this case, the peak at 
4662.0 //Hz would be due to noise. 

The echelle diagram with the fourteen identified modes is 
shown in Fig. ^ The frequencies of the modes are given in 
Table |3 As the curvature of the p-mode alignment could be 
important above 5.5 mHz, both identified high-frequency £=1 
modes can take two different values shifted by 1 1 .57 //Hz. This 
should be taken into account when comparing individual fre- 
quencies with theoretical models. In the same way, the peak at 
5475.7//Hz was identified as a ^ = 2 mode, taking a probable 
curvature at high frequency into account, but could just as eas- 
ily be a radial mode. The large separations are given in Fig.|5] 
and have a mean value of 161.7 + 0.3 //Hz. 

3.3. Oscillation amplitudes 

Concerning the amplitudes of the modes, theoretical compu- 
tations predict oscillation amplitudes near 13 cms"' for a star 
like 70 Oph A, with mode lifetimes on the order of a few days 
(Houdek et al. 1999 1. The amplitudes of the highest peaks are in 
the range 1 l-14cm s"', in agreement with the theoretical value 
and with the observations of a similar star, a Cen B (Carrier & 
BourbanEnOD- 

4. Conclusion 

Our observations of 70 Oph A yield a clear detection of p-mode 
oscillations. Several identifiable modes appear in the power 
spectrum between 3 and 6 mHz with an average large spac- 
ing of 161.7//Hz and a maximal amplitude of 14.5 cms"'. The 
mono-site observations, coupled to the low resolution of the 
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Table 3. Oscillation frequencies (in /iHz) for the large spacing 
of 161.7yuHz. The frequency resolution of the time series is 
2.2 //Hz, and n' is the radial order defined within a constant 
value. Another possible identification can be done by inverting 
all I - 0,2 and (- \ modes, except for the peak at 4662.0^(Hz, 
which would be due to noise in this case. 



n' 


C = Q 


e= 1 


1 = 2 


-7 


3378.1 






-6 


3540.5 






-5 








-4 


3862.1 






-3 








-2 


4182.8 


4258.8 




-1 


4343.7 















1 






4662.0 


2 








3 


4993.2 


5066.0 




4 


5153.7 






5 
6 


5316.9 




5475.7 


7 
8 




5717.8 
5879.3 





4000 6000 6000 

Frequency [^tillz] 

Fig. 5. Large spacing Av versus frequency for p-modes of de- 
gree t-Q (•) and (-1 (a). Open symbols correspond to large 
spacing averages taken between non-successive modes. All in- 
dividual errors are fixed to V2 x l.lyuHz (half resolution). 
The dashed-line corresponds to the mean large spacing of 
161.7yuHz. 
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time series and to the faint signal-to-noise, do not allow us to 
unambiguously disentangle {-2 from {-Q modes. The small 
spacing seems to be consistent with theoretical predictions for 
such a star with a value near lO^Yiz. However, the lack of cou- 
ples of ^' = and 2 modes for the same radial order does not 
give a single identification solution: £ - 0,2 can be replaced by 
{=1 modes and ^ = 1 by ^ = 2 modes. The study of the 70 Oph 
system, with asteroseismic and non-asteroseismic constraints, 
has been postponed to a second paper. 
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